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The chromoelectric polarizability of mesons governs the strength of the gluonic van der Waals force 
and therefore of non-quark-exchange processes in hadronic physics. We compute the polarizability of 
heavy mesons with the aid of lattice gauge theory and the Born-Oppenheimer adiabatic expansion. 
We find that the operator product expansion breaks down at surprisingly large quarks masses due 
to nonperturbative gluodynamics and that previous conclusions concerning J/T/i-nuclear matter 
interactions and J/ip dissociation in the quark-gluon plasma must be substantially modified. 



I. INTRODUCTION 

Although hadronic interactions are a central phe- 
nomenon of nuclear and hadronic physics, very little can 
be said about them from first principles and a micro- 
scopic description remains elusive. One of the few at- 
tempts at describing hadronic interactions from QCD 
originated more than twenty years ago with the oper- 
ator product expansion (OPE) approach of Peskin 1]. 
Peksin argued that the coupling of soft external glu- 
ons to small (heavy quark) hadrons may be considered 
a short distance phenomenon and therefore amenable 
to Wilson's operator product formalism. More recently, 
Luke, Manohar, and Savage have placed Peskin's ar- 
gument in the context of effective Lagrangians. Briefly, 
in the absence of flavour exchange processes, the inter- 
actions of hadronic matter are dominated by multigluon 
exchange processes which may be described with an effec- 
tive Lagrangian at the compositeness scale, Ag ^ fq^ ^ 
as{AQ)mQ: 

^eff - E -^iPivYPiv)-V,{vyVivr){cEOs+CBOB). 

V Q 

(1) 

Here the gluonic operators are Oe = —G^^^'G^v^vp and 
Ob = iG'^^G^p - G^'"G^^v^vp■, P{v) and Vf,{v) create 
pseudoscalar or vector mesons with four velocity v'^. In 
the meson rest frame these operators reduce to E'^ and 
respectively. The multipole formalism becomes exact 
in the large quark mass/small hadron limit and the cou- 
pling constants (or Wilson coefficients) ce and cb may 
be interpreted as the chromoelectric and magnetic polar- 
izabilities of the heavy meson. 

Luke, Manohar, and Savage used this formalism to es- 
timate the binding energy of quarkonium (T and ^') with 
nuclear matter The gluonic matrix elements were es- 
timated with the aid of the scale anomaly and the exper- 
imentally determined gluonic momentum fraction of the 
nucleon at the scale Ag. The final ingredient was Peskin's 
original estimate of the chromoelectric polarizability 
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(in the large Nc limit). We note that theoretical uncer- 
tainty is introduced through the choice of the compos- 
iteness scale, the strong coupling ag, and the size of the 
meson. The final estimates of the binding energies were 
roughly 3 MeV for the T and 10 MeV for the J/-0. Subse- 
quently, Brodsky and Miller 3J used this result to obtain 
a J/i/'-nuclear matter scattering length of as = —0.24 fm 
and a cross section of roughly 7 mb at threshold. Brod- 
sky and Miller also argued that multiple gluon exchange 
dominates the J/\E'-nuclear matter interaction. 

Finally, Kharzeev and SatzQ have applied Peskin's re- 
sults to the interaction of 7/4" with comoving matter in 
heavy ion collisions. They argue that the cross section is 
small near threshold and that therefore collision-induced 
dissociation should not confound the use of J/vP suppres- 
sion as a diagnostic for the formation of the quark-gluon 
plasma. 

The chromoelectric polarizability has appeared in at 
least one other context. Leutwyler has argued that non- 
perturbative level shifts in the heavy quarkonium spec- 
trum may be related to the product of the vacuum expec- 
tation value of the electric field pair density and the elec- 
tric polarizability!^ (similar arguments have been made 
with QCD sum rulesjal)- In particular he states that the 
small size of the heavy meson implies that quarks inter- 
act with slowly varying random chromofields. The energy 
shift is then given by the expectation of the operator 



SH = -PE • r- 



1 



Ha — Es 



-E-rP 



(3) 



where P projects onto mesonic states which are orthog- 
onal to the meson, E^ is the mass of the heavy meson, 
E is the chromoelectric field, and Ha is the Hamiltonian 
which describes the interactions of quarks in the colour 
octet state, 



Ha = 2mQ 



as 



2N.r 



(4) 



CE 



3(iV2 - 1) 



(2) 



The potential in Ha is the perturbative expression for the 
interaction of a quark and an anti-quark in the colour 
adjoint representation. The expectation value of 6H is 
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proportional to ce{E'^) and thus the electric polarizabil- 
ity gives the strength of nonpcrturbative mass shifts (or, 
equivalently, the strength of the nonlocal nonpcrturba- 
tive potential due to interactions with the gluon conden- 
sate). 

It is clear that the value of the electric polarizability 
is crucial to all these conclusions. In the following, we 
carefully examine Peskin's computation of ce and con- 
clude that its true value is roughly a factor of ten smaller 
than claimed. More importantly, it will be shown that, in 
this application, the operator product expansion is never 
reliable in Nature due to the effects of nonpcrturbative 
gluodynamics. 



II. CHROMOELECTRIC POLARIZABILITY 

Peskin specifies two conditions which permit the ap- 
plication of the operator product expansion. The first is 
that the meson should be small, t'q^ 3> ^qcDi which im- 
plies toq KqcdI o.s[rQ^). The second constraint arises 
because gluons coupling to the heavy meson must arrange 
themselves into colour singlets. Thus the emission of a 
single gluon - which raises the energy of the meson to 
that of an octet (or hybrid meson) state - must be fol- 
lowed quickly by a subsequent emission. The correlation 
time between these events is At ~ l/(£'a — E^) where Ea 
is the energy of the intermediate hybrid state. Thus the 
colour singlet criterion is Ea — E^ ^ es ^ ^QC d which 
imposes the stronger constraint: 



mQ-> KqcdI al- (5) 

We have introduced the Coulombic binding energy ts — 
mQClallA {Cf = {N^ - l)/(2iVc)) and the principle 
quantum number of the heavy meson, n. Eq. |S1 implies 
that the potential felt by the heavy quarks is perturbative 
and hence that the heavy meson wavefunction is nearly 
Coulombic. Peskin estimates that the condition of Eq. O 
is met for mg ^ 25 GeV for n ~ 1. We note that this 
result is obtained in the large Nc limit where Ea tends to 
2mg and hence At ^ l/es. An updated limit in which 
this constraint is considerably relaxed will be established 
in section IV. 

Under the conditions specified above, gluon emissions 
must arrange themselves into small colour singlet clus- 
ters which are attached to a small region in spacetime in 
which the heavy meson is in an octet state. This obser- 
vation permits the application of the operator product 
expansion. Peskin applies this idea by exponentiating all 
possible two-gluon couplings to the heavy meson. The re- 
sult is a gauge invariant effective interaction of the form 

T _ Sr^ JN)i] 3 2-2N T^iT^2N-2T^] (f.\ 

^eff - - 2^(^E o-o<^B ■ E h-> (b) 

where Dq is the temporal component of the covariant 
derivative. We follow Peskin and introduce the dimen- 



sionful parameters oq — 21 [CFdsi^Q) (the Coulombic 
Bohr radius) and es to make the Wilson coefficients, 
dimensionless. We note that the leading term has 
already been given in covariant form in Eq. ^ and as a 
model in Eq. 13 

The expression for the Wilson coefficient is 
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and (j) represents the heavy meson of interest. For S- 

wave states c^'*'"' = 5'^^c^^\ Finally, using Is Coulombic 
wavefunctions and neglecting the adjoint potential yields 

the result0 for c^^(ls) given in Eq. [21 (we suppress the 
superscript from now on). A similar computation gives 
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c_b(1s). 



(8) 



In general rg ^ Kruqas), and the energy denominator 
scales as mga^/n^, thus CE{ns) ^ nP (factors of mg and 
as cancel against the prefactors in Eq. [Tj) . It is clear that 
the OPE breaks down very quickly with the principle 
quantum number. 

As we have already remarked, these estimates have 
been used to compute the strength of hadronic interac- 
tions in a variety of applications. However, a number 
of strong assumptions have been made in deriving them. 
Certainly, it is not clear that the adjoint potential need 
be as simple as perturbation theory indicates. Fortu- 
nately recent improvements in lattice gauge theory have 
allowed for an accurate determination of this interaction 
in the heavy quark regime [3- It is therefore expedient 
to confront the assumptions of Refs. Q, II, H, |^ IE Q with 
lattice gauge theory in an attempt to establish the valid- 
ity of Eqs. [5] and Thus we briefly review the current 
knowledge of hybrid potentials before moving on to a re- 
evaluation of the polarizability and the operator product 
expansion itself. 



III. ADIABATIC HYBRID SPECTRUM 

A simple consequence of the fact that glue is confined 
is that it must manifest itself as a discrete spectrum in 
the presence of a static colour source and sink. In this 
case the physical hadrons are heavy hybrid mesons. It 
is relatively easy to study heavy hybrids by constructing 
gluonic configurations on the lattice which are analogous 
to those of a diatomic molecule. Indeed, the gluonic con- 
figurations may be described with the same set of quan- 
tum numbers as diatomic molecules: A^. Here the pro- 
jection of the total gluonic angular momentum onto the 
QQ axis is denoted by A which may take on values E, H, 
A — 0,1,2, etc. The combined operation of charge and 
parity conjugation on the gluonic degrees of freedom is 
denoted hy rj = u, g and Y ~ ± represents reflection of 
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the system in a plane containing the QQ axis. As with 
the diatomic molecule, all systems with A greater than 
zero are doubly degenerate in Y . Gluonic adiabatic sur- 
faces may be traced by allowing the heavy quark source 
and sink separation to vary and hybrid mesons contain- 
ing excited gluonic configurations may be studied in the 
adiabatic Born-Oppenheimer approximation. 

The results of a recent lattice computation are pre- 
sented in Figure The lowest state is the S+ surface 
and corresponds to the Wilson loop static interquark po- 
tential. The first (second) excited state is the n„ (Ilg) 
surface and may be visualized as a gluonic flux tube with 
the addition of a single 'phonon'. A similar analogy exists 
for all of the higher states. 
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scales. Thus the constraint is not nearly as strong at 
finite Nc and with reasonable hybrid potentials. One 
concludes that gluons are largely correlated in time as 
required, lending hope to the idea that the application of 
the OPE to hadronic interactions may be unexpectedly 
robust. Unfortunately, a new constraint exists, which we 
now demonstrate. 

A central criterion for the validity of the operator prod- 
uct expansion is that the hybrid surfaces shown in Fig- 
ure ^ approach a universal form at short distances. Fig- 
ure n makes it clear that a universal hybrid potential 
behaviour (namely Vg) does not appear until 

K(rQ) >yA^'(rQ)-FA^(rQ). 

Since the typical hybrid surface separation is order Aqc d 
for small tq (and is much larger for the splitting relevant 
to ground state mesons), one has ttiq 3> ^h.QCD 1 0?^ ~ 
150 GeV. Thus Eq. [Sjis recovered (albeit with an un- 
lucky additional large factor); however, the constraint is 
now a necessary condition rather than merely sufficient 
as before. Thus, although the condition which insures 
the emission of colour singlet states is likely to be satis- 
fied for all quark masses, the hidden assumption in the 
method, namely that a universal octet potential is rele- 
vant, is only true for very heavy quarks. 

This conclusion has a simple interpretation in 
Leutwyler's random field model: the appearance of a 
discrete hybrid spectrum makes it clear that the correct 
representation (here we employ the Born-Oppenheimer 
approximation) of the matrix element of SH of Eq. Ois 
as follows 



FIG. 1: Low Lying Adiabatic Hybrid Surfaces. Lattice data 
for the (circles), Ilti (diamonds), and Hg (squares) surfaces 
from Ref. Q. Lines are simple parameterizations of the data. 
The scale is ro ~ 1/2 fm. Inset; Potentials at Short Distance. 
The same data is shown in traditional units. The IIu and 
Ilg lines are those of the main figure with the addition of the 
perturbative adjoint potential Va- 

The inset of the figure shows the lowest hybrid surfaces 
at scales less than 1.5 femtometres. The dashed and dot- 
ted lines are those of the main figure with the addition 
of the perturbative adjoint potential, Va — Q!s/(6rQ) joj]. 
The figure demonstrates that perturbative behaviour has 
not been seen at rg « 0.2 fm (points) and that this is 
consistent with expectations (dotted and dashed lines). 



SEn = {c^n\SH\(t>n) 

^ (0„;S]+|(5i/|0„;E+) 

^ |(0„;I]+|E-r|fe;A^)P 

In this expression h represents all of the nongluonic quan- 
tum numbers which describe an intermediate heavy hy- 
brid state. The essence of the operator product expan- 
sion is that this expression factorizes (ze., the Wilson 
coefficients depend on short range physics only). Factor- 
ization requires that the hybrid energies in the denomi- 
nator do not depend on the gluonic quantum numbers, 
A^. It is only in this circumstance that the expression 
simplifies: 



IV. DISCUSSION 

We first note that the nonperturbative gluodynamics 
shown in Figure ^ indicates that the octet-singlet split- 
ting used in Eq. [S] is not accurately described by the 
Coulombic binding energy. Rather, the figure indicates 
that this splitting is roughly 1 GeV at typical hadronic 



'^r.-Ejl^yi^t\^'K) (10) 

and the operator product formalism is valid. 

Finally, we consider the value of the chromoelectric po- 
larizability in light of the lattice hybrid data of Figure Q] 
We choose to numerically compute the (/> wavefunction 
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in the Born-Oppenheimer approximation with the aid of 
the lattice surface. The sum over intermediate hy- 
brid states is performed numericaUy by expanding in the 
eigenstates of the Ilg surface (this is the lowest surface 
which couples to a vector heavy hybrid - which is the 
case we consider in the following). 

The results are shown as the open squares in Figure |21 
One sees that the Peskin result of Eq. |21is recovered fax- 
row) in the very heavy quark mass limit, as expected|lOj. 
However, the value of C£;(ls) at the T or J/^P masses 
(arrows on the abscissa) is highly suppressed with re- 
spect to the asymptotic value. The diamonds are numer- 
ically obtained values for the polarizability in the case 
that the adjoint potential has been included in the Ilg 
surface parameterization. Again, the analytical result, 
C£;(ls) = 2347r/425'55, is approached very slowly in quark 
mass. 




FIG. 2: Chromoelectric Polarizability as a Function of Quark 
Mass in GeV. Points represent ce(1s) as computed with the 
Ilg surface with (diamonds) and without (open squares) the 
perturbative adjoint potential, Va- The line is the approxima- 
tion of Eq. ^3 Inset: The Ratio ce{2s) /ce{1s) versus Quark 
Mass. 



It is tempting to speculate that the majority of the 
finite quark mass correction to Peskin's result is due to 
the hybrid mass gap. Allowing for this in Eq. [7| yields 
the following generalization of Eq. |21 



Stt 



— 256(1 + u)^/^ -256 



3(iV2 - 1) 
384f - 96v^ + 16v^ - 6v^ + 3v' 



(11) 



where v 



and Vb is the strength of the rel- 
The resulting 



evant hybrid potential at its minimum, 
expression is shown as a solid line in Figure |21 evidently 
the agreement is quite good and this expression may serve 
as a useful extrapolation to light quark masses. 

Finally we display the ratio C£;(2s)/ce(1s) in the inset 
of Fig. 121 Again the Peskin result (Eq. |SJ) is approached 



only very slowly in the heavy quark limit. We note that 
Eq. |S1 leads to the uncomfortable prediction that the T' 
interacts 5000 times more strongly with nuclear matter 
than does the T. The inset shows, however, that this pre- 
diction is substantially moderated (from 5000 to roughly 
15) when finite quark mass effects are taken into account. 



CONCLUSIONS 



According to the arguments of Refs pj,|j,|5|,|6j, the elec- 
tric polarizability of a small meson controls the strength 
of its interactions with hadronic matter via the operator 
product expansion. We have recomputed this strength 
with the aid of lattice hybrid potentials and find that 
the large mass gap between the ground state (S+) and 
excited state gluonic configurations leads to a strong 
suppression of the electric polarizability as the quark 
mass is reduced. The result is that, if one neglects is- 
sues of the applicability of the OPE, previous estimates 
of interaction strengths are reduced by roughly a fac- 
tor of 100. Thus the arguments of Khazeev and Satz 
concerning the utility of J/ip suppression as a quark- 
gluon plasma diagnostic are strengthened. Alternatively 
non-quark-exchange J/^/^-nuclear matter interactions are 
greatly reduced, suggesting that quark exchange mecha- 
nisms should be carefully considered in the analysis of 
the J/T/i-nuclear matter binding issue. 

On a more general level we have argued that the as- 
sumptions underlying the operator product expansion de- 
scription of heavy hadron interactions are violated for all 
physical states. This situation arises because the suffi- 
cient condition on temporal correlations among giuons 
has been replaced with a necessary condition on the ap- 
plicability of factorization which is only true for very 
heavy quarks. 

The short length scale required for factorization arises 
for a number of reasons. Certainly the fact that the 
strong coupling is not large and that the ratio of funda- 
mental and adjoint Casimirs is also small, help to under- 
mine the reliability of the OPE. However, the relatively 
flat behaviour of the adiabatic hybrid surfaces below 1 fm 
must be considered the leading cause. One may speculate 
that this arises due to the robust persistence of string- 
like field configurations, even at quite small interquark 
separations. It thus appears that strong nonperturbative 
gluodynamics conspires to bring about the demise of the 
operator product formalism in this application. 

Although we have said nothing about the utility of 
the OPE in the very heavy quark limit, the authors of 
Ref. 8] show that the interaction between very small 
colour dipoles becomes nonperturbative (it is essentially 
correlated two pion exchange). It thus appears that the 
premise of the OPE and any effective field theoretic ap- 
proach to the interactions of small hadrons is compro- 
mised. 
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